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Abstract 



X 

c^ I Given a space it is easy to obtain the system of geodesic equations on it. In this paper the 

inverse problem of reconstructing the space from the geodesic equations is addressed. A 
procedure is developed for obtaining the metric tensor from the Christoffel symbols. The 
procedure is extended for determining if a second order quadratically semi-linear system 
can be expressed as a system of geodesic equations, provided it has terms only quadratic 
in the first derivative apart from the second derivative term. A computer code has been 
developed for dealing with larger systems of geodesic equations. 



1. Introduction 



The development of geometry originates in its applications for map making [1], but even 
more from its use in kinematics and dynamics [2]. As such, it is of interest to look at 
the interplay of geometry and dynamics. The path of a test particle in a flat space is 
a straight line. The curved space generalization of the straight line is a geodesic. Thus 
test particles in curved spaces move along geodesies. Given a space one easily obtains the 
system of geodesic equations on it [3]. In principle it should be easy to obtain the space 
from the geodesic equations. This would be of physical relevance as actual observations 
would not provide the space but would provide the paths followed by "test particles". 
For example, in general relativity, one assumes some matter-energy distribution and then 
solves the Einstein equations to obtain the metric tensor, using which one obtains the 
geodesic equations, which give the paths for test particles [4]. However, one does not 
really know the matter-energy distribution in any actual situation, but only the observed 
paths of particles. Consequently, it would be of interest to be able to determine the metric 
directly from the geodesic equations. 

Though simple in principle, the problem is complicated by the fact that the ChristofFel 
symbols are non-linear combinations of the metric tensor and its first derivative. As such, 
a system of highly non-linear first order partial differential equations would have to be 
solved to obtain the metric tensor. The problem can be reduced enormously in complexity 
by contracting the Christoffel symbols with the metric tensor, to obtain a system of first 
order linear partial differential equations. In general, even this is very complicated to 
solve. Further, we would need to check compatibility of the solutions obtained. 

The procedure adopted here uses the skew symmetry of the covariant form of the Riemann 
tensor in the first two indices, and the symmetry under interchange of the first and second 
pair of indices to provide a system of linear equations that can be solved simultaneously. 
If the system does not decouple we finally have to solve n partial first order differential 
equations for one function of n variables. Consequently an arbitrary constant appears in 
the solution. If it decouples we need to solve correspondingly more differential equations 
and hence more arbitrary constants appear. These constants generally get determined by 
inserting the solutions back into the equations for the metric. In principle it is possible 
that they may not be fully evaluated and lead to a class of metrics. In this paper a 
specific prescription is provided to reconstruct the space from the geodesic equations. 
The metric constructed is unique (up to some multiplicative constants appearing in the 
solution). For two variables the procedure will be explained in detail in the next section. 
However, for larger systems the procedure is still too complicated to be implemented by 
hand. In section 3 we have provided a general discussion of the general case and the 
logic of the computer code to obtain the metric from the Christoffel symbols. Further, 
if we only have a system of second order quadratically semi-linear ordinary differential 
equations (ODEs) given, that have only the quadratic term, we would not know whether 
they could, consistently, be regarded as a system of geodesic equations. In section 4 we 
give a brief discussion of how the code can check whether the system can, or cannot, be 
regarded as a system of geodesic equations. In section 5 there are some specific examples 



given to illustrate the use of the general procedure. In the last section we have given a 
brief summary and discussion of the results. 



2. A System of Two Equations for Two Variables 



The essential principle for obtaining the metric from the Christoffel symbols may be 
seen directly by considering a system of two geodesic equations for two functions of one 
variable. However, the general procedure involves additional complications that will be 
discussed later. For the system of two equations 

x" = a{x, y)x''^ + 2b{x, y)x'y' + c{x, y)y''^, (1) 

y" = d{x, y)x''^ + 2e(x, y)x'y' + /(x, y)y''^, (2) 

we can read off the Christoffel symbols as the the negative of the coefficients of the 
quadratic terms. Thus 

J- 11 = ~0'-, J- 12 ~ ~") -'- 22 ~ ^^1 -'- 11 ~ ~"'i -'- 12 ~ ~^i -*- 22 ~ ^/- \^) 

Note that for a general second-order quadratically semi-linear system of ODEs, the coef- 
ficients cannot be assumed to be expressible as Christoffel symbols. However, if we are 
given the system of equations as geodesic equations, we can assume that the coefficients 
are so expressible. For a known metric tensor the Christoffel symbols are then given by 

r}fc ^ ■^5'*"'(fi'jm,fc + gkm,j — gjk,m)- (4) 

Now construct the Riemann tensor from these Christoffel symbols 

r)i -pi -pi I -pi -pm -pi -pm /r\ 

jkl ~ ■*■ jl,k -*- jk,l "T -■- ml^ jk ^ mk^ jl 5 I'JJ 

where the Einstein summation convention, that repeated indices are summed over, has 
been used. Note that the tensor is skew in the last two indices, k and /. As such, when 
they are equal the tensor is trivially zero. Thus, without loss of generality, we can set 
k = 1,1 = 2. Putting the tensor into fully covariant form it is skew in its first two indices 
as well. Using the metric tensor to lower the index of the curvature tensor, we obtain the 
two linear relations for the metric coefficients: 

^ll^\l2 + 9l2R\l2 = 0, (6) 

fi'12-R 212 + fl'22-R 212 = 0- (7) 

There are various possibilities for the Riemann tensor components being zero or non-zero. 
Not all possibilities are consistently allowed. Apart from the case of a flat space, -R*^^ = 0, 



only two possibilities survive: (a) when i = j, BJ'jti = 0; (b) when i = j, R^jm 7^ 0. In 
case (b) these equations can be used to write gu and (722 in terms of gi2 

R^ 

911 = —Br^9i2 ■■= Agi2, (8) 
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^22 = -^^12 := Bgu. (9) 
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To make the procedure easier to see, write gu = p{x, y),gi2 = (l{x, y) and (722 = ^'{x, y). In 
case (b) using eqs.(7) and (8) we get p and r in terms of q. Then, writing the Christoffel 
symbols explicitly we obtain the differential equation for q 

q^ = {Ab + a + Bd + e)q, qy = {Ac + b + Be + f)q. (10) 

The solution for q is provided by integrating eq. (9) relative to x and y and comparing the 
arbitrary functions of integration 

q{x,y) = a{y)exp{ {Ab + a + Bd + e)dx) = (3{x)exp{ {Ac + b + Be + f)dy) . (11) 

There would appear to be an arbitrary constant still left. This disappears on using the 
resulting p, q, r in the expression for the Christoffel symbols. (Remember that the inverse 
metric contains the functions as well as their first derivatives.) 

In case (a) g = 0. We now get two sets of two partial differential equations for p and r, 
which can be solved to give 

p{x, y) = '^{y)exp{ / 2a(x, y)dx) = 6{x)exp{ / 2c(x, y)dy); (12) 

r{x, y) = fi{y)exp{ / 2d{x, y)dx) = u{x)exp{ / 2/(x, y)dy). (13) 



There now appear to be two arbitrary constants appearing which are determined by 
inserting the expressions back into the Christoffel symbols. If the constant (s) remain, we 
obtain a class of metrics for the same geodesies. 

Note the remarkable fact that the algebraic symmetry properties of the 'geometric' in- 
trinsic curvature tensor, from the 'ODE' point of view, are just the compatibility criteria 
for being able to obtain the metric from the system of geodesic equations. 

Finally, there remains the case that the space is fiat. The metric tensor can certainly 
be set as the fiat space metric tensor in Cartesian coordinates and hence the metric is 
"reconstructed". However, this would not be the metric tensor in the coordinates used. 
We could now solve the full system of six linear first order partial differential equations 
for the three functions p, g, r of two variables x, y. The compatibility is now guaranteed. 
There are other, neater, methods available as well [5]. 



3. The General Procedure 



Whereas, in principle there is nothing new when we have more than two variables, the 
problem arises because there are many more possibilities now. To see how the equations 
proliferate, consider the three variable case. We now have 18 Christoffel symbols for a 
system of three geodesic equations. These lead to six independent components of Rijki- 
There are now three sets of constraining equations each of which has three possible index 
choices. As such, we have 9 linearly dependent equations for the 6 metric coefficients. If 
all the components are non-zero, we have enough equations to be able to obtain the metric 
coefficients from three differential equations. In fact if we have five distinct components 
non-zero, we could solve the system. However, we have many possibilities between this 
case and the flat metric. For 4 variables there are 10 independent metric coefficients, 
20 linearly independent components of Rijki and 40 Christoffel symbols for a system of 
four geodesic equations. This appears to be a very heavily over-determined system and 
compatibility checks would become really long. 

The proliferation of equations would have rapidly rendered the problem intractable were it 
not for the availability of computer codes to solve such systems, for many more variables, 
going through all possibilities. We have constructed such a computer code that enables 
us to solve the problem in full generality. It is given at: www.cam.wits.ac.za/inverse. 

The logic of the code is as follows. We first differentiate the Christoffel symbols r*j^ 
relative to all the dependent variables and combine them to form the curvature tensor 
of valence [1, 3], i.e. with one upper and three lower indices, R^jki- Next we use the 
symmetry properties of the covariant form of the Riemann tensor, namely 

QimR jkl = ^9jrnRikli (14) 

and 

QimR jkl = QkniRlij- (15) 

Since (14) is skew in i, j, there are n^(n — l)/2 linearly independent equations. Further, 
(15) are n^ equations. There are only n{n+ l)/2 independent components oi Qij. As such, 
the system must be grossly over-determined. However, if the F*^ are, indeed, Christoffel 
symbols, they must be consistent. As such, one can use the first n{n+ l)/2 — 1 of them to 
obtain all the Qij in terms of one of them (say gn). It is to be noted that since the system 
is homogeneous, there can be no non-trivial determination for all the metric coefficients 
from here. One now writes the equation for the metric tensor in terms of the given 
Christoffel symbols as 

dikj + gjk,i — gij,k = '^dii^jk- (16) 

With the full set of equations for all 'n?{n+ l)/2 independent Christoffel symbols we can 
reduce the equations to a system of n first order linear partial differential equations for 
one function (say gn) of n variables. We can now solve these and obtain the full metric 
tensor. 



It may happen that the system of equations has a rank less than n(n + 1)/2 — 1. If it is one 
less, we will need to solve the system for n{n + l)/2 — 2 components of the metric tensor 
and then solve partial differential equations for the last two metric coefficients. Similarly, 
if the rank is p less, we would solve the system for n{n + l)/2 — p — 1 and then solve the 
remaining partial differential equations for those p components. Note that there would be 
no need to re-check compatibility of the solutions, other than to determine the arbitrary 
constants arising from the solution of the differential equations. 



4. Consistency Criteria for Systems of Geodesic Equa- 
tions 



So far we have taken it for granted that the system given is for geodesies. It may happen 
that one obtains a system of equations that look formally like the system of geodesic 
equations, in that they can be written as 

x' + Ti^x^x'' = 0, (17) 

but that they cannot be regarded as a system of geodesic equations. The point is that there 
is no check that the system of partial differential equations (16) is internally consistent. 
To check this we would require that (14) and (15) form a consistent set of equations. This 
is still not enough). We also need to check that the first Bianchi identities are satisfied, 
namely 

R'jki + R\ij + R\jk = 0- (18) 

If these are satisfied we can, indeed, regard the given system of equations as a system of 
geodesic equations and possibly use the results of theorems on global linearizability of the 
system to obtain the solution [5]. 

The computer code we have prepared can be used to obtain the metric tensor if the 
system is known to be of geodesic equations and can be used to check the consistency of 
the system as geodesic equations. 



5. Computation 

The algorithm is implemented as follows: Specify the order of the geodesic equation via 
n. Assume Qij = Qji and -Rj„„ = by enforcing the following rules 

SetAttributes[g, Orderless]; (19) 

R[L, j_, k_, L] := 0; i == j&&k == 1; (20) 



in Mathematica. We next introduce the lists SkewSymmetry and Symmetry which 
caters for all the possible combinations of i,j,k and I when summing over the repeated 
index in (fl^ and (TTSl) . The independent variables are represented by X in eql[i_, j_, k_, L] 
and eq2[i_, j_, k_, L]. The matrices which are then formed when mapping eql[i_, j_, k_, L] 
and eq2[i_, j_, k_, L] to their corresponding lists Symmetry and SkewSymmetry , are 
stored in SYM and SkewSYM respectively. We construct the gij metric tensors with 
gcomponents , which is then used in conjunction with CoefRcient Arrays to construct 
the matrices ASym and ASkewSym. By choosing n{n + l)/2 — 1 rows from each using 
the input from ChooseEqns these two matrices are used to form the matrix A and vector 
b. The vector sol then solves the metric tensors in terms of gn by using LinearSolve in 
conjunction with the matrix A and vector b. The overdetermined system of linear partial 
differential equations, EqnSetlG , are then used to solve for gu by using DSoIve. The 
order of the problem dictates the number of arbitrary functions which will then have to 
be solved subsequently. 

The implementation of the n = 3 case is given in the Appendix. 

The code has been checked for the following examples. 

1. Systems of two equations: (a) geodesies on a sphere; (b) a linearizable system 
given in [5]; (c) a non-geodesic system (in [5]). 

2. Systems of three equations: (a) flat space; (b) a 3-sphere; (c) linearizable (in [5]); 
(d) non-geodesic (in [5]). 

3. Physical four dimensional Lorentzian systems: (a) the Reissner-Nordstrom 
system (geodesies for a point charged mass, in which the charge could be taken to be zero 
to get the Schwarzschild system); (b) the Kerr system (for a rotating point mass, in which 
the rotation could be taken to be zero to reduce to the Schwarzschild system). 



6. Summary and Discussion 



We have shown explicitly how to construct the metric from the geodesic equations. In 
other words, if we knew the geodesies globally we could reconstruct the full manifold 
with the metric on it and if we know them locally we can reconstruct the metric and 
hence the space, locally. It is remarkable that the purely geometric entity measuring 
curvature should provide, when looked at from the viewpoint of differential equations, 
the compatibility conditions for the system to be regarded as describing geodesies. The 
significance of this representation is that it provides a procedure to reduce the system 
of equations as follows. We use the conjecture of [6] that for m-dimensional sections 
of constant curvature it will have an so{m + 1) symmetry algebra. We can use the 
geometric information to choose the sections of constant curvature to decouple the system 
of geodesic equations (using the procedure of [6]). These m geodesic equations would then 
be completely solved and we would have only a system oi n — m coupled equations to be 
solved. Notice the heavy utilization of purely geometric considerations. 



We solved the problem explicitly for the case of a system of two geodesic equations ana- 
lytically. The only problem arose in the case of the flat space where we anyhow know the 
metric in Cartesian coordinates. If we want to express the metric in the given coordinates, 
it may not be so easy. However, the metric coefficients are directly determinable by solv- 
ing decoupled first order partial differential equations for each of the metric coefficients. 
There are also more elegant methods available [5]. 

For larger systems one needs a computer code. Such a code was developed and its logic 
is given here. The code is available at: www.cam.wits.ac.za/inverse. Some examples 
illustrate the use of the code. The code further provides a check for the given system of 
equations to be consistently regarded as a system of geodesic equations. 

This approach is of importance as it provides a geometric method for solving systems 
of ODEs [5]. Further developments may be possible by embedding the space in higher 
dimensional spaces in which they become larger quadratically semi-linear systems. At- 
tempts to convert cubically semi-linear systems to quadratically semi-linear systems in a 
higher dimension are in progress [7]. This procedure is the inverse of that adopted by 
Aminova and Aminov [8], in which they use the geodetic re-parametrization symmetry 
(d/ds) to reduce the system by one dimension. Further, it would be possible to reduce 
the order of a higher order system by increasing the number of variables. Thus we could 
possibly use the same techniques for higher order ODEs by embedding in correspondingly 
higher dimensions. This technique could also be tried to reduce from higher degree equa- 
tions to two. It would be of great interest if the approach could be extended to PDEs as 
well. 



Appendix 

We implement our code for n = 3. The skew symmetry of i?*;,; implies that without loss 
of generality, in both (Ell) and ([221) below we have {k,l) = {(1,2), (1,3), (2,3)}, while in 
([22D (i,j) = {(1,2), (1,3), (2, 3)}. The equations (^ and dSSj) are our skew symmetry 
and symmetry equations repsectively and are generalized as 

9jiRiki + 9iiRjki + 9j2Riki + 9i2Rjki + 9j?.Riki + 9iiRjki = 0' (21) 

9ilRjkl - 9klRuj + 9i2Rjkl - 9k2Rlij + 9i3Rjkl - 9k3Rlij = 0. (22) 

From ( 12T|) we obtain the nine equations which is generated by eq2[i_, j_, k_, L] 

5'11-Rll2 + 9l2Ru2 + 9l?,Rll2 = 0; (23) 

giiRl,, + guRlis + 9i3Rlu = 0, (24) 

911RI23 + 912RI23 + 913RI23 = 0, (25) 

guRlu + ^22i?2i2 + 923RI12 = 0, (26) 



S'12-R213 + 922R2I?, + fi'23-R213 = 0' (27) 

5'12-R223 + 922R22Z + fi'23-R223 = 0' (28) 

5l3i?312 + ^23^312 + 9^iRll2 = 0, (29) 

5'13-R313 + ^'23-^313 + ^'33-^313 = 0' (30) 

^'13-^323 + ^'23-^323 + fl'33-R323 = 0- (31) 

From ( |22|) we obtain the six equations which is generated by eql[i_, j_, k_, L] 

5'ii (^213 - ^312) + 9i2 [Rli3 - Rlu) + 9i3 [Rli3 - RI12) = 0) (32) 

9llR22'i + 9l2 (-R223 ~ -R312 j ~ 922R312 + fi'l3-R223 ~ fl'23-R312 = 0) (33) 

9ii (-R312 - -R213) + 9i2 (-R312 - -R213) + fi'13 (-R312 - -R213) = 0) (34) 

9llRz2Z + 9l2 (-R323 - -R313) ~ fi'22-R3i3 - 5^23-^313 + 5^13-^323 = 0' (35) 

911R22Z ~ 9i2 [R312 ^ -^223) ~ 922R312 + 9nR223 ~ 92ZR312 = 0' (36) 

- ^'11-^323 + 9l2 (-R313 - -R323) + ^'22-^313 + ^'23-^313 ^ ^'13-^323 = 0- (37) 

We choose fl2^ . fl27|) . fl?I]) . ([32D and ([33]) to solve for for gi2, gi3, g22, 923 and g^^. This 
choice is done by ChooseEqns = {4, 5, 9, 10, 11} in our code, which is then used to 
relate these metrics to the gn metric. Here are the relations 

^ ~ -"^212-^^223 \R2i3) ~ 

pi ( p3 \ I pi p3 p3 I ( p2 p3 I p2 p3 \ p3 
"-'^212 y-'^SUJ ' -"'212-'^213-'^312 ' \^-'l'312-'l'212 "•" -'l'212-'l'213y -^^223 

O Q/Q q\\o "1 0/'q\ 

+it223-fL2i2 (^-"^213 ~ -"^312 jj -"^213 ~ -'^213-'^212 \R3i2) + 

p3 I p2 p2 p3 I p2 f pi p3 I pi p3 _i_ p2 p3 | A 1 

-'T'213 \^-n'212-'T'223-'^213 ">" "'^312 \^ -'T'213-"'212 "r -"-2 12 -"-2 13 "•" -'l'212-'l'223y j "•" 

-r^312 (-^^213-^212 ^ R212R213) \R2i3 ~ R312) + 

(-^213-^^312-^212 + (-^^212 \R2i3 ~ R223) " -"^212-^^312 j -"^213J -^312' (38) 

5'12 = 5'11 ((-^213-^212 - -R212-R213) ((-^312 " -^213) -^223 + -^223 (-^213 ~ -R312))) /^^ (39) 

fi'13 = fi'll ((-R312-R212 - -^212-^312] -^213 + 

p2 p2 I pi p2 \ p3 I p2 / p2 p3 _i_ pi p3 \ _i_ 
''^9.T?.-'^29.3 ' -'''91 9 -'''3 19 -'''913 ' -'•'913 I -'•'993 -'•'91 9 ' -'•'9 1 9 -'•'3 1 9 ' 



212-"'223 ~ -"'212-"'312y -"'213 ~ -"'213 \^ -"'223 -"'212 ~ -"'212-"'312 
-R3I2 (-R9I9 (-R91?i + -Rsi9 ~ (-R9IS + -Rsi9 -R9I9 -^913 + 



312 V-fl'212 ^-^^213 "T -"'312^ ~ ^-'^'213 "^ -'^312^ ^212 j j -"'213 

p2 f p2 p3 p2 p3 \ I 

-"'312 \^ -"-213 -"-212 -"'212-"'213y "•" 

9 



-^223 (-^213 ~ -R312 j (-^212-^213 ~ -^213-^212) + 
-R312 ( (-^212-^223 + -R212-R312 R2 



1-212-' 1-223 ~ -'i-212-'i'3l2y -"1-213 
-^213 (-^223-^212 + -R212-R312)) ) /^) (40) 

fi'22 = fi'll (^213-^212 - -R212-R213) ((-^213 " -R312) ^223 + ^223 (-^312 " -R213)) /^) (41) 

fi'23 = -fi'll (-R213-R212 - -R212-R213) ((-^213 - -R312) -^223 + -^223 (-^312 " -^213)) /^) (42) 
^33 = fi'll ((-^212-^323-^223 + -^323 (-^212-^312 ~ -^312-^212)) (-^213) + 

(^-"^323 (-n^212-n^213-"^223 + -"^223 -"^212 (-^^213 ^ -"^312 jj + 

-^323 ((-^212-^223 + -^^212-^^312 j -^213 ~ -^213 (-^223-^212 + -^212-^312 jj + 

-R312 (-R323 ((-^213 + -R312J -^212 ^ -^212 (-^213 + -^312J j ~ -^212-^323-^223)) -^213 + 

(-^312) -^323 (-^212-^^213 '~ -^213-^212 j + -^223 (-^323 (-^213 "~ -^312 j (-^213-^212 



p2 p3 \ I pi p2 p2 / p3 p3 

-rt2i2-fi-213i "'"-'^-212-^^213-^323 1-^1-213 -^^312 



pi / £?! d2 p2 p3 I pi / p2 I d2 p3 _i_ pi p3 \ 

-'1-312 \^-'i-212-'^-213-'l-323-'l-223 "•" -'^323 \^-'l-213 \^-n'223-'l-212 ">" -'l-212-n'312 j 

,^212^223 + ^212-^312) ^213)) ) /^- (43) 



Substituting into (16) we obtain a system of differential equations for the i?k; in terms 
of the Christoffel symbols which is represented in our code by EqnSetlG. Imposing 
Christoffel symbols we can determine the i?k; from the system of ordinary diferential 
equations and hence the metric. On solving for gu we get 

fill (a:i,X2,X3) = (44) 

fill (a:i,a;2, 0:3) = 



(^S2 {K2, X2, Xs) Tl^ {K2, X2, X3) + S3 {K2, X2, X3) Vf^ {K2, X2, X3)) rfi^2 

We can then subsequently solve for the arbitrary function Ci [x2,a;3]. 
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